Swapping entangled Kondo resonances in parallel-coupled double quantum dots 
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Strong electron and spin correlations are studied in parallel-coupled double quantum dots with 
interdot spin superexchange J. In the Kondo regime with degenerate dot energy levels, a coherent 
transport occurs at zero temperature, where two entangled (bonding and antibonding) resonances 
are formed near the Fermi energy. When increasing J or the dot-lead parallel-coupling asymmetry 
ratio r2/ri, a swap between two entangled resonances occurs and the line shapes of the linear 
conductance are interchanged. The zero-bias differential conductance shows a peak at the critical 
values. Such a peculiar effect with the virtue of many-body coherence may be useful in future 
quantum computing. 



A large number of proposals have been made to ma- 
terialize quantum bits (qubits) and quantum computing. 
Among these proposals, coupled quantum dot (QD) sys- 
tems are particularly attractive. 1] The spin degree of 
freedom of the localized electrons on the dots is con- 
sidered as a qubit due to the comparatively long co- 
herence time. A key challenge is the construction of 
coupled double QD (DQD) to perform a swap opera- 
tion, i.e. exchanging the electron spin states on the 
two dots. When the square root of a swap operation is 
combined with other isolated qubit rotations, a quantum 
controUed-NOT gate can be built and any quantum algo- 
rithms can be implemented. ^2] In the coupled DQD, two 
local electrons form a singlet state. It has been proposed 
that the swap operation can be realized by tuning the 
time-dependent interdot spin superexchange (ISS) J(t) 
from positive to negative, flopping the singlet and triplet 
states. Ill 



In this Letter, we propose a simple and reliable mecha- 
nism to perform such a swap process in a parallel-coupled 
DQD at low temperatures. It has been well-established 
that under the Coulomb blockade with odd number elec- 
trons on a single QD, a quantum coherent many-body 
(Kondo) resonance is formed near the Fermi energy in 
the dot density of states (DOS) 0|. The Kondo effect for 
even number electrons on a single multilevel dot and pos- 
sible phase transitions between singlet and triplet states 
have been considered for both "vertical" [sj and "lateral" 
Q, S, configurations. For the coupled DQD with 
degenerate energy levels, the Kondo behavior and the 
ISS interplay and strongly cornpete, as seen in the bulk 
two- impurity Kondo problem, [lOj and a question arises 
whether there exists a spin entangled state composed of 
the coherent Kondo resonances. '11 Jr^revious experimen- 
tal and theoretical studies have ma inly focused on the 
serial-coupled DQD 0,El[ll|ll|ll[llHII3: ex- 
cept for Ref. 0,0. However, it has been recently realized 
that the serial geometry is unsuitable for studying this 
competition experimentally and a direct evidence for ob- 
serving spin entanglement between the dot local electrons 



could be sought in the parallel coupled configuration. ..2fl] 
It has thus motivated us to investigate whether and how 
this competition manifests itself in the coherent trans- 
port through DQD in the parallel configuration. 

For a degenerate DQD with ISS, a special dot-lead 
coupling configuration (see Fig. ^ is considered. By in- 
creasing the dot-lead coupling asymmetry ratio r2/ri 
(ri.2 = '^PfVi2)y where is the DOS at the Fermi level 
and Vi,2 is the dot-lead hopping integrals, the dot local 
electrons couple to the right and left electrodes can be 
transformed from serial (r2 = 0) to symmetric parallel- 
coupled configuration (r2 = Fi). We generalize the slave- 
boson mean field (MF) theory [ilj to take into account 
both electron and spin correlations simultaneously. Us- 
ing an entanglement order parameter (EOF) Af and the 
asymmetry ratio r2/Fi to describe the entanglement be- 
tween the local electron spins, we find that two entangled 
bonding and antibonding Kondo resonances are formed 
very close to the Fermi energy. When increasing the ISS 
or the dot-lead coupling asymmetry ratio, the EOF can 
change sign at a critical value. As a result, the bond- 
ing and antibonding resonances swap, or equivalently, 
the singlet and triplet levels interchange and the zero- 
bias differential conductance displays a peak at the crit- 
ical value. Importantly such a swap effect occurs only 
when the dot-lead parallel couplings are asymmetric, i. e. 
< r2/Fi < 1 and the gate voltage controlling the in- 
terdot electron hopping is fine tuned. 

We describe the coupled DQD system as two Ander- 
son magnetic impurities with degenerate levels and in- 
finite on-site Coulomb repulsion, and an antifcrromag- 
netic(AF) spin superexchange between two local electron 
spins is generated by the second order perturbation in 
the interdot electron tunneling. The model Hamiltonian 
is given by 

k,CT;a=L,_R o-;i=l,2 
k,cr 
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FIG. 1: Configuration of the coupled double quantum dots 
with asymmetric couplings to the left and right electrodes. 
Two degenerate dot levels (ei ~ £2) axe coupled by an AF 
spin superexchange with strong on-site Coulomb repulsion. 



(Va/iV^i + Vifl^b^) C^,a,R + h.c. 
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where N is the total number of electrons, the slave- 
boson representations di^a = bifida have been used to 
describe the local electrons on each QD {bi and fi^a de- 
note the respective hole and electron occupied states), 
and the local constraints bjbi + J2a fiafi,cr = 1 have to 
be imposed j^]. 81 and 82 correspond to the spin den- 
sity operators of the dot local electrons, characterized 
by = E...' 4.^^,.' with r« (a = x,y,z) the 
Pauli matrices. Here, we are only concerned with the 
equilibrium (ek.L = ek,fl) properties. A similar coupling 
configuration of the DQD model has been studied in the 
nonintcracting case. 22] Since the ISS interaction can be 
written in an SU(2) singlet form 0| 

cr,o"' a, a' 

in the MF approximation, an EOF Aj — J2aifi af^-cr) 
can be introduced to describe the spin singlet between 
the local electrons on two separate dots. As in the usual 
treatment of the Anderson impurity model at T = 0, 
the bosonic operators are replaced by their expectation 
values and the local constraints by the Lagrangian mul- 
tipliers Xi. It has been established that the slave-boson 
MF treatment captures the basic Kondo physics for the 
single-impurity Anderson model at low temperatures, 
and such a theory becomes exact for large local spin 
degeneracy, [i^ In the presence of degenerate dot energy 
levels, we expect to have Ai = A2 = A and bi — b2 = bo, 
namely, the degeneracy of the dot energy levels can only 
be lifted by the ISS interaction. Thus, an effective model 
Hamiltonian is obtained 



k,cr;a 
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E [/^^ (^k,.,L - Ck,.,/i) + h.. 



where e^d = + A and Vi = b^Vi are two renormalized 
parameters, and two canonical modes of bonding and an- 
tibonding have been introduced by ao- = (/i,cr + f-2,a) 
and (}^ — ^ {fi,a ~ /2,(t)- The difference between 
the bonding and antibonding energies are just given by 
IJAf = Jeff- Actually, Jeff also corresponds to the 
effective energy splitting of the spin singlet and triplet 
states formed by the renormalized local electrons on the 
dots. Thus, a singlet-triplet transition can occur when 
the EOF Af changes sign, instead of reversing the sign 
of J. 

When a Nambu spinor <i>J. = (/i,<t,/2,ct) is defined, 
the Fourier transform of the retarded Green's function 
-(TT-$^(r)$t (r')) can be derived as 



uj ~ ed + «(ri + F2) 




JAf + 2i\lTiT2 















where Gx is the Pauli matrix, and the DOS on each 
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QD is given by Ariijj) 

corresponding to a superposition of the bonding and 
antibonding resonances lying at energies Iq,/? — (id T 
JAf) with renormalized hybridization widths Ta,f3 = 
bgTrpj: {Vi ± ¥2)^- Actually, these are two entangled reso- 
nances with many-particle coherence, and we will refer to 
them as entangled bonding and antibonding Kondo reso- 
nances. At the symmetric parallel-coupling r2 = Pi, only 
the bonding combination of the conduction electrons, i.e. 
(Ck,CT,L + Ck,CT,fl) couples to the localized electrons, and 
the antibonding combination (Ck.o-.L — Ck.o-.fl) are com- 
pletely dropped out. 

To determine the MF order parameters fog and A / , the 
corresponding self-consistent equations '^„{fi„fi,(y) — 
l-bl and J2a{fi,af2,cr) = A/, are rewritten as 



1-bl 




(2) 
(3) 



where the first equation represents the total quasipar- 
ticle occupations on bonding and antibonding energy 
levels, while A^ corresponds to the difference between 
these two occupations as follows from the Friedel sum 
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Moreover, a third self-consistent equation from 
= {[bi,H]) =0 is needed, yielding 



(a) 



Fa In 



D 



Tp In 



D 



(4) 



where two characteristic energy scales are formally de- 

fined as r„(^) - ^ (e^ T JA/)' + {hl^Pf? {Vi ± ^2)^ 
with D as the bandwidth of the conduction electrons in 
leads. We have to stress that these two energy scales 
can not be directly related to the Kondo temperatures 
defined in the multilevel QD system. To find the 
saddle point solution numerically, we choose Fi = 1 as 
the energy unit and D — 100. 

In the Coulomb blockade regime for each dot, the de- 
generate dot energy level is set by = —6, far below 
the Fermi energy. We find that grows as increasing 
while T[3 is a small constant (T^ < T^). For a given 
value of both Ta and Tp are almost independent of 
J. In the parallel-coupled DQD with a tunable ISS, the 
electron occupation on each dot strongly depends on pa- 
rameters F2 and J. Given a dot-lead coupling asymmetry 
ratio r2, we find that the MF value of A/ grows as J in- 
creasing and changes sign at a critical value Jc- When 
A f reverses sign, the relative position of the bonding and 
antibonding energies of the renormalized dot energy lev- 
els are switched, i.e., the quantum spin states on QDs 
are swapped. Similarly, for a fixed value of J, there also 
exists a critical value of the dot-lead coupling asymme- 
try ratio (0 < r2,c < 1), where A/ changes sign. Such 
a sign change effect signals the existence of two different 
regimes in the asymmetric parallel-coupled DQD system. 
In Fig[21 we display EOF A/ as functions of the ISS J 
and the ratio of the dot-lead couplings r2. Actually, the 
dramatic effect of Ay is not sensitive to the parameter 
as long as the condition \€d\ > T2 is satisfied. In Fig[2t, 
a phase diagram is given in the three-parameter space 
(ed, r2, J). Moreover, we evaluated the dot electron DOS 
for ed = —6 and F2 = 0.6 with different ISS, shown in 
the left cohmm of Fig|3| When J = 0.1, A/ has a very 
small value, and the entangled bonding and antibonding 
resonances are very close to each other and hardly distin- 
guishable in Fig|3^. For J = 0.205, A/ reaches its largest 
value and the entangled bonding resonance appears be- 
low the Fermi energy (FigEt). When J = 0.210, A/ 
changes sign and the bonding resonance state swaps its 
position from below to above the Fermi level. At J = 0.3, 
the bonding and antibonding resonances emerge again. 
In this process, a sharp and narrow entangled resonance 
always stands close to the Fermi level. 

From the Landaucr formula, the linear conductance 
can also be calculated 



2e^ 



Tr [G/(w - iO+)rflG/(w -I- iO+)ri] , (5) 
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FIG. 2: The order parameter A/ in the Kondo regime (ed = 
—6) as functions of the asymmetry ratio of the dot-lead cou- 
plings r2 in (a) and the ISS J in (b). The ground state phase 
diagram is given in (c). Fi = 1 is as the energy unit. 



where G/(cj + iO+) and G/(tj - iO+) 
sent the retarded and advanced Green' s 



tions, and Tj 



= 2fe2 
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func- 
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are two matrices re- 



lated to the asymmetric dot-lead couplings. The 
corresponding results are delineated in the right column 
of Figl^l For the intermediate couphngs of J, the fincar 
conductance displays a Lorentzian conductance peak 
centered at the bonding energy and a Fano resonance 
at the antibonding energy (FiglSJ^'). The latter arises 
due to the presence of bound states of DQD embedded 
in the conduction band continuum. However, when 
the EOF Af changes sign, these two characteristic line 
shapes are switched (seen in Fig|3t'). For a smaller 
value of J in Fig|3K', the Fano resonance with a small 
transmission is an anti-resonance at the Fermi energy, 
because of the destructive quantum interference between 
different pathways through QDs. As J is large in 
FigEli', there is a progressive increase of the width of 
the bonding resonance. For the serial-coupled (r2 — 0) 
and symmetric parallel-coupled (r2 = Fi) DQD, these 
swapping features do not appear. To make contacts 
with experiments directly, we extract the zero-bias 
differential conductance G{u! = 0) as functions of the 
asymmetric coupling parameters F2 and ISS J for a 
given value of — —6, displayed in Fig^Jt and 0Jd. 
Surprisingly, we find that the differential conductances 
have maxima precisely at the critical couplings for swap. 
The appearance of such sharp conductance peaks can 
be explained by the majority occupation of the bonding 
quasiparticle states across the Fermi level during the 
swap process. We should point out that the singlet- 
triplet transition discussed here is very different from 
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FIG. 3: The dot DOS and the corresponding hnear conduc- 
tance for different interdot AF spin superexchange J with a 
given asymmetry ratio of the dot-lead couplings r2 = 0.6 and 
degenerate dot energy levels = —6. 
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FIG. 4: The zero-bias differential conductance as functions of 
the asymmetry ratio of coupling parameters r2 and the ISS 
J with degenerate dot energy level ed = —6. 



what was considered in multilevel dots |a, LDj ISl l2| j where 
the coupling parameter symmetry excludes completely 
the antibonding channel of the conduction electrons, so 
the swap effect does not appear there. It is true that the 
quantum fluctuations beyond the MF description may 
change detailed behavior near the swap point, but the 
existence of the swap itself is a robust effect, because 
two MF solutions with opposite sign of A/ are stable. 

To conclude, in a degenerate parallel-coupled DQD 
with asymmetric parallel couplings, two entangled bond- 
ing and antibonding resonances are formed close to the 
Fermi energy in the Kondo regime of each dot. A swap 
effect between two resonances has been found, leading 
to a sharp peak centered at the critical coupling in the 
zero-bias differential conductance. To observe such a pe- 
culiar effect in experiments, one has to fine tune the cor- 
responding gate voltage controlling the interdot electron 



hopping or the dot-lead coupling asymmetry ratio in a 
parallel-coupled DQD system. Moreover, this effect will 
lead to a practical and reliable mechanism to construct 
the quantum gate for the future quantum computing. 
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